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Abstract. Under appropriate assumptions the energy of wave equations with 
damping and variable coefficients c(x)utt — div(fe(a;)VM) + a{x)ut = h(x) has 
been shown to decay. Determining the rate of decay for the higher order 
energies involving the fcth order spatial and time derivatives has been an open 
problem with the exception of some sparse results obtained for fc = 1, 2, 3. We 
establish estimates that optimally relate the higher order energies with the 
first order energy by carefully analyzing the effects of linear damping. The 
results concern weighted (in time) and also pointwise (in time) energy decay 
estimates. We also obtain L°° estimates for the solution u. As an application 
we compute explicit decay rates for all energies which involve the dimension 
n and the bounds for the coefficients a(x) and b{x) in the case c{x) = 1 and 
h{x) = 0. 



1. Introduction 
We will study hyperbolic equations of the form 

(1.1) c{x)uu-Ai\{h{x)Vu)+a{x)ut^h{x,t), a; £ R", t > 0, 

where a,b,c, and h satisfy the assumptions listed below. This system is generally 
accepted as a model of wave propagation in a heterogeneous medium with friction 
(given by a{x)ut) and source terms h{x^ t). The coefficient c{x) accounts for variable 
mass density, while h{x) is responsible for temperature changes as the wave travels 
in space (see the derivation in [6]). More surprisingly, (jLlj) has been considered as 
a model for heat propagation by Cattaneo [2] , and independently by Vernotte [17j . 
In this setting the constitutive equation for the flux q as given by 

(1.2) qt+a{x)q^b{x)yu 

replaces the classical Fourier's law q = b{x)Vu. The Cattaneo- Vernotte equation 
(|1.2p is usually considered as appropriate to describe unsteady heat conduction. For 
more on the applicability of this model see [31 [HI [S] ■ The finite speed of propagation 
feature that the solutions enjoy under this formulation renders this system as a more 
accurate way to study the real effects of heat diffusion. We will show later how these 
aspects regarding the physical interpretation of (|1.1|) tie in with the mathematical 
analysis of finding the decay rates, but we believe that this paper offers compelling 
evidence as to why the diffusion phenomenon is an essential feature for the damped 
wave equations. 
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The literature is replete with results on energy estimates for hyperbolic systems 
with linear damping, but they mainly concern the case when b{x) — c{x) = 1 and 
h{x,t) = 0. In [13j we considered wave equations with variable coefficients in the 
case c{x) — 1, h(x,t) = for which we obtained explicit decay rates that depend 
on the dimension and the growth of the coefhcients. The main goal of this work is 
to show how the energies are related to each other; more precisely, we deduce decay 
rates for the higher order energies in terms of the decay rate for the first energy 
E{t). This issue was addressed for example in [TT] where the author shows that for 
c{x) = b{x) — 1, and h{x,t) — the energy has the following decay rate 

Eit;u) :=^i\\u,it)\\^ + \\Vu{t)r)<{l+tr\ 

while for the second order energy the decay rate improves so one has 

E{t;ut) = liWuttim' + \\^u,{t)\f) < il + t)-\ 

By using these estimates in the equation one consequently obtains 

\\Au{t)\\'<{i+tr\ 

The idea of using the diffusion effect in order to show a faster decay for the higher 
order energies is also present in ^ where the authors show that for the elasticity 
system with linear damping ut one has 

E{t; u) + E{t; Ut) + E{t; Uu) < (1 + t)-\ 

One expects, however, that the gain in the decay rate when moving from the 
first energy to the second should be t~^. This is motivated by the fact that the 
long time behavior of the linearly damped wave equation with constant coefficients 

(1.3) Utt — Au + auf — 
resembles the behavior of the corresponding parabolic equation 

(1.4) - Au + aut = 0, 

as it was suggested in [10], [12], [16]. By estimating the higher order energies for 
the Gaussian as a solution of (|1.4p . one notices that E(t; d^^^u) decays faster than 
E{t; d^u) by a factor of t^'^. This suggests that a similar phenomena should occur 
for (jl.Sp . Indeed, our paper shows that this conjecture is correct, so in the absence 
of other effects caused by inhomogeneities, one basically has 

E{t; d^u) < E{t, u){l + t)-2fe, A: = 0, 1, 2, ... 

In this spirit, by matching the gain in decay t-^ of E{t] dt+\) versus E{t; d^u) 
obtained for diffusion problems we show the optimality of our decay results for 
higher order energies. One of the most important contributions of this paper is 
developing a methodology for determining the rates of decay for energies of all 
orders for a variety of hyperbolic systems with linear damping. 

The paper is organized as follows. In the next section we present the assumptions 
for our results, after which we present the main results regarding the decay of the 
higher order energies involving time derivatives. The weighted decay of the spatial 
derivatives is obtained for orders fc = 2 and fc = 4 from which we deduce L°° 
estimates of the solution in dimension n = 3. We apply these results in the case 
c = l,/i = in section 6. In the Appendix one can find the theorem and proof 
regarding the finite speed of propagation of solutions of (|l.ip for variable h and c. 
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2. Main Assumptions 
Let a, 6, and c be smooth coefficients satisfying the conditions: 
ao(l + |x|)-" <a(a;) <ai(l + |a;|)-", 

(2.1) boil + \x\f <b{x) <biil + \x\f, 
co{l + \x\)-^ <c{x)<ci{l + \x\)-\ 

where a^, 6^, q, i = 0, 1 are positive constants and the exponents a, 7 satisfy: 

(2.2) 2-/3-7<2, 2a + /3-7<2. 
We consider initial data 

(2.3) u\t=o^uoeH\R"), ut\t=o=ui€L^{R''), 

and a source h e L°°{[0, 00), L^(R")). The three fmictions have compact support: 

(2.4) uo{x)^0, ui{x)=0 for \x\ > R, 

(2.5) h{x,t)^0 for q{\x\) > t + q{R) , 

where q which measures the increase in the support of u is given in the Appendix. 



3. Weighted L'^ estimates for time derivatives 

The strategy behind determining the rate of decay for E{t; d^u) in terms of the 
decay for E{t; u) is based on some simple facts. First, since d^u solves the equation 
(j2.3p with the different inhomogeneity d^h, the higher order energy E{t;dfU) sat- 
isfies (almost) the same decay estimate as E{t; u) with possibly different constants 
(the actual decay will depend on d^h). The question is whether these decay rates 
can be improved similarly to the case a = b = c = 1 and h = where the Fourier 
representation of u shows faster decay rates of norms involving higher derivatives. 

Higher order spatial norms are estimated from 

div(6(a;)Vu) = a{x)ut + c{x)utt — h(x,t), 

since taking directly cc-derivatives produces new terms and changes the dissipative 
nature of this equation. In other words, we will express the x-derivatives in terms 
of the ^-derivatives. 

Let us recall the definition 

E{t; m) = i y" {cul + 6| Vwp) dx 

and energy identity 

(3.1) —E{t;u)-\- [ aui dx — [ hutdx. 

dt J J 

A simple consequence is the following. 

Proposition 3.1. Assume that < To < T and > 0. Then 

I I {l + tfa{x)u'^tdxdt < {l + TofE{To;u)+ f [I + tY'^ E{t;u)dt 

JTq J JTq 



To 



,„h^ix,t) , , 
(l + tY \ \ ' dxdt. 
a[x) 
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Proof. The energy identity (|3.1|) implies 
^[(1 + t)f'E{t; u)] + (1 + t)^ / au^^dx = (1 + t)^ I hutdx + + t)''-'^E{t; u). 



dt 

Notice that \hut\ < + /2. The resuh foUows from integration on [ro,T]. 

□ 

We can now proceed with higher-order norms. Let v — ut io simphfy notations. 
Then 

c{x)vtt - div{b{x)'\/v) + a{x)vt = /it(x, i), x G R", t > 0, 

and 



(3.2) 



-— j {cvl + b\V v\'^) dx + J avj dx = J htVtdx. 



In order to obtain integrals of E{t; v), we multiply the equation for v with ■^W{t)v 
where W is to be chosen later. Integration on R" yields 

1 d f Wa-Wtc 2\ , f „r 2 , 
-— [WcvtvA — —v ] dx- Wcv^dx 

(3.3) +\ j [Wcv^, + Wb\yu\^ + ^»^-^^% 2^ rf^ 

\ I WhfV dx. 



2 

We will choose W satisfying 

(3.4) W{t)< inf Wttc-Wta>0. 

xe supp u{-,t) c[x) 

The support of the solution u at time t is given by 

suppu(-,i) = {x G M" : |x| < R + Ct^^} 

for some C > 0, as a consequence of the Corollarv lA.2l from the Appendix. 
Adding and we have 

-— / {cvf + b\Vv\^ + Wcvtv-\ — -V j dx+ {a-Wc)v^ dx 

(3.5) +i I (Wcv^ + Wb\Vv\' + WuC~W,a ^^^ 
htVt + ^Whtv^ dx. 

Assume that W also satisfies 

(3.6) < ciif + b\Vv\^ + Wcvtv + ~ < 2cvf + b\Vv\^ + CoWav'^, 
where Co is a constant; later we give explicit conditions in terms of W., a and c. 
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Next multiply identity (|3.5p by {1 + 1)" and integrate on [TbjT]: 



To 



+ / [cvt + b\\'v\^ + Wcvtv- 



Wa - Wtc 



dxdt 



1 

2 

f-T 



To 



(1 + tf {Wcv^ + Wb\Vv\^) dxdt 



< 



Hence, 



'To 



(1 + t)" ( htVt + -Whtv ] dxdt. 



(1 + ty {Wcvf + VF&I Vwp) dxdt 



Ut + 5|Vt;p + Wcvtv + 



Wa - Wf C 



dx 



t=To 



2 



To 



(1 + t)''-' cvi + felVi-r + WcvtV 



Wa-Wtc 2 



V dxdt 



/To 



(1 + ty ( htVt + -WhtV ) dxdt. 



Our assumptions on the quadratic form yield a simpler estimate: 

r-T 
/To 



(1 + ty {Wcvt + Wb\Vvf) dxdt 



< {i + nr 



{cvl + 6|Vvp + CaWav'^) dx 



2\ 



+1^ / (1 + i)''"^ (ctj^ + b\Vv\^ + CoWav^) dxdt 



{1 + ty {2htVt + Whtv) dxdt. 



'To 



Choose To sufficiently large, such that 
(3.7) {1 + tyW >2iy{l + ty-^ fort>To, 

so by combining similar terms we obtain 



1 



'To 



(1 + ty {Wcv'i + Wb\\/v\^) dxdt 



(3.8) <{l + Tay 

I To 



{cvf + b\Vv\^ + CoWav^) dx 

■ T 

To 



t=To 



+iyCo I f{l + ty-^Wav^ dxdt + f f{l + ty {2htVt + Whtv) dxdt. 
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We can estimate the last integral as follows: 

' (1 + tf {2htvt + Whtv) dxdt 



(3.9) 



< 



+ tycWv^ + C(e)(l + ty{cW)-^h^t] dxdt 

dxdt 



i(l + ty-^Wav^ + i(l + ty+^Wa-^h^t 



where < e < 1/4. 

Choose a smooth function W such that it satisfies 

(3.10) wi{l+t)-'^ <W{t) <w2il+t)-'^ 

for some < a; < 1 and wi,W2 > 0. This implies that 

(1 + tYicwy^h^t < (1 + ty+^wa-^hl 

which together with (|3.9p simplifies (|3.8p to the following: 



(1 + ty {Wcvl + Wb\\/v\^) dxdt 



(3.11) 



<(l+To) 



[cvf + fclVup + Wav'^) dx 



t=To 



+ 



To 



[l + ty^^Wav^ dxdt+ / I {l + ty+'^Wa-'^hi dxdt. 

To 



The growth assumptions from (j3.10p used in this last inequality yield: 



To 



{l + ty^'^ {cv^+b\Wvf) dxdt 



(3.12) <(i + ro)'- 

cT 



{cv'^ + 6|Vup + Wav'^) dx 

■ T 



"-"^-^av^ dxdt 



To 



t=To 

(1+0 



^^uj+l^-lf^2 ^^^^ 



By using the estimate from Proposition 13.11 with fi = v — lu — 1 and (|3.4p we 
derive the following estimate: 

(3.13) f I {1 + ty-"^ E{t; v) dxdt < {1+ Toy [E{To;v)+E{To;u)] 

J To J 

+ [ {l+ty-'^-^E{t;u)dt+ [ [{l + ty-'^-^a-\{l+tfh^t+h^)dxdt. 



To 



To 



The proof is complete provided we show: 

Existence of the weight function W. Let 

(3.14) W{t) ^wo{l+ty 

where : 

(i) wi < wq < W2 with wi, W2 given by (|3.10p . 
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(ii) the exponent ui is chosen such that 
(3.15) maxio, ~ i> < ^ < 1. 



2-/3-7. 

This choice for the weight W satisfies all the constraints as we show below: 

• The inequality p.4p i . By (|2.ip it suffices to show: 

Wa{l + t)-^<— inf ^ (l + |x|)^-". 

\x\<R+Cf^~l^-i 

For a < 7 this is obvious since (1 + — > as t ^ oo, while in the right 
hand side we have (1 + > 1. For a > 7 we need to show 

Cl 

but for sufficiently large times t this holds since lj > ^'^^^ by (I3.15p . 

• ([01) 9 holds since for W given by (|3.14p we have Wtt > 0, Wt < and the 
coefficients a and c are positive. 

• For the left inequality in (|3.6p we first complete the square in v and vt, so 
we are left to show 

2Wa - 2WtC - W'^c > 0. 

By using ()2.1|) and ()3.14|) we reduce the problem of proving the above 
inequality to showing: 

(3.16) 2woao{l + ty^{l + |a;|)-" + 2woUJCn{l + ty'^-\l + \x\y 

-«;gci(l+t)-2-(l + |a;|)-T >0. 

This can be simplified to 

(3.17) 2w;oao(l + |x|)-"+'^ + 2u;owco(l + t)-^ - woCi{l + t)-'^ > 0. 

The assumption p.lSp on lu gives us that the above inequality is true for 
large times. 

We follow a similar approach for the right inequality in p.6p ; we complete 
the square so we have: 

v,~—j +Co—v - —v>0. 

We need to show 

(3.18) {4:Co-2)Wa + 2WtC-W^c>0 

The constant Cq can be chosen sufficiently large so that we have 4Co — 2 > 0. 
By (|2.ip and (|2.2p in order to prove (|3.18p it is enough to show 

(4Co - 2)woao{l + \x\)-"+'' - ciwl{l + ty^ - 2wwoCi(l + t)-^ > 0, 

which is done exactly as above when proving p.l7p . 

• The condition (|3.7p is satisfied since w < 1. 



Denote hy 9 ^ v — uj. From p.l3p with a simple induction argument we obtain 
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Proposition 3.2. Let a, b, c, and h be smooth coefficients which satisfy S2.1]} and 
i2.^) . and let 9 > 0. Then the solution u of equation with initial conditions 

\2. Sl\) . which satisfy \2.4-^ , satisfies the following weighted energy estimate: 

(3.19) / {l + ty+^''E{t-d1u)dt<{l + nY[Y,E{n;dlu)\ 



'^0 i=0 



+ f {l+tfE{t;u)dt+f I {l + tY+^a-^Y}^^ + ^)'''^^thf]dxdt. 

The above arguments allow us to also obtain pointwise decay rates for the energy 
as stated in the following: 

Proposition 3.3. Under the assumptions of Proposition [37^ we have 

k 

(3.20) E{T;dS) <{1+Ty 



{l+TonY,EiTo;diu)] 



+ /^l + tfEit; u)dt + f /(I + tr^ E ^l±&^d.dt 

"'To "/To J a\X) 

Proof. In (|3.ip written for we apply Young's inequality with appropriate coef- 
ficients such that the damping term disappears. This yields 

(3.21) E{T- d'^u) < E{t- d'^u) + J < * < ^■ 
Hence 

(3.22) / + t)^+^''E{t; d^u)dt > E{T; d^u) / (1 + tf+^'^dt 

"^'(1+t)-- r [ \^^^d.dsdt. 

To Jt J aix) 

Simplifying the above inequality and using (|3.19p we obtain the desired estimate. 

□ 

Let us go back to lower order derivatives and estimate 

aufdx = J —cuttUt + div{bVu)ut — hufdx 

which by the Cauchy inequality gives 

/ r \ 1/-^ / /■ \ 1/2 

,2 ^ I „„,2 i„ 1 I „ ,2, 



(3.23) J aUf_dx < cu^^dx j u^dx 

(3.24) + (^J b\Vu\'^dx^ (^J b\Vut\^dx^ + ^ J '^dx + ^ J au^^dx. 
We obtain the following 

Proposition 3.4. Under the assumptions of Proposition [37^ we have 

(3.25) / au^dx < {E{t;u)E{t;ut)y/^ + / —dx. 
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The explicit decay can be obtained by using (I3.20p with fc = and k — 1; the 
decay rate of energy associated with the damping term is approximately —1 over 
the decay rate of total energy, in agreement with Proposition 13. II on the "average" 
in time decay rates. 



4. Weighted estimates for spatial derivatives 

Decay estimates of higher order spatial derivatives are important for studying 
nonlinear perturbations of equation Such problems require more regular 

solutions and naturally lead to W'^'^ norms with fc > and 2 < p < oo. Although 
the latter norms are expected to decay faster than H'^ norms, this is not easy 
to verify. The first difficulty is that norms are not convenient to estimate 
by the multiplier method when p ^ 2. Here the standard approach is to use 
interpolation or embedding, i.e., the Gagliardo-Nirenberg or Sobolev inequalities. 
We thus reduce the question to multiplier estimates with losses of derivatives. The 
variable coefhcients present an additional difficulty as x-differentiation produces 
terms that change the dissipative form of p.f p . A simple solution is to express 
x-derivatives in terms of ^-derivatives from the equation: 

c h 

Mu = -utt + ut , 

a a 

where Mu — a^^div(6VM). The diffusion phenomenon means Mu w Ut, so second- 
order spatial derivatives are related with the first-order time derivative. Similarly 
we can derive an identity for M'^u. 

Below we give two sufficient conditions on a, 6, and c to guarantee decay estimates 
of Mu and M^u, respectively. Given Ai, A2 G [0, 1], the coefficients satisfy 

2 

Proposition 4.1. Let Mu — a^^ div{bVu) and assume that u is a solution of 
(i) If 173) holds, then 

a{Muf < (1 + t)^^cu^^ + auj + —. 

a 

(a) If (gJp and 114. 2\ ) hold, then 

a{M^uf < (f + tf^^cul^, + (1 + tf^ {cul, + b\Vuu?) 
+ {l + tf-aul 

+M + (i + <)2AiM* +afAf- 
a a \ a 

Proof, (i) The first claim follows from (|4.f p and 

/ /j2 

aiMu)"^ < 3a —u^^ + u^ + — 



(4.1) 
(4.2) 



sup 

xG. supp u{- .t) 



sup 

xG supp u(-,t) 



^ b{x) 
a{x) a{x) 

1 



,c(x) 



a{x) 



div b{x)\7 



a{x) 
a{x) 
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(m) To verify the second claim, we apply M{uv) = uMv + vMu + 2(6/a) Vm ■ Vi). 
We have the chain of identities 



f c h 

M -utt + ut 

\a a ^ 

uttM- + -Mutt + 2-V- • Wutt + Mut - M-. 
a a a a a 



Using the expression for Mu, we further obtain 

, c c 



UttM- + - ( -utttt + uttt - — 
a a \ a a 



+2- V- • Vutt + -Uttt + utt ] -M- 



a a 

9 

c c 
= -^utttt + 2-Uttt 



a 



a 



V a 



+ l]utt + 2-V- -Vutt 
a a 



a a a a 



The square of M^u is bounded by the sum of squares times 7: 



a{M'^uf < 7a 



—Utttt + 4— w 



a 



2 "ttt 



M- + 1] u 

a 



'la 



iVuttI 



2 , h\ 



c2 hi 



M- 



Finally, we rewrite last estimate to match conditions (|4.1|) and (|4.2|) : 



a{M'^uf 

Thus, we have 

a{M'^uf 



< 



a 



V- 

a 



b\Vutt\' + ^ + ^-^+a{M- 



< 



(1 + tf^'cul,, + (1 + t)^i (c^.?,, + b\Vutt\^) 



hi 



+ (1 + ty^auit + — + (1 + tf^^ — + a ( M- 



a 



□ 



It is possible to derive expressions for M'^u, fc > 3, and find conditions on a, 
b, and c which yield decay estimates of higher order spatial derivatives. However, 
considering fc < 2 is sufficient for most applications in M" when n < 3. 

Proposition 4.2. Let a,b,c, and h he sufficiently smooth functions which satisfy 
( fO) and Define 

divibVu) 
Mu = i '-. 
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Then the following weighted estimates hold for the solution u of equation fHiP with 
initial conditions 112. 3\) which satisfy (2^: 



('0 



{l + ty+'aiMuy dxdt 



'To 



< 



1 X 

il + Tof+'[J2EiTo;dlu)]+ [ [l + tf E{t-u) dt 



To 



{l + tY\dlh) 



{l + tf+^-^^a{M'^uf dxdt 



I To 



< 



3 T 

il + Tof+'[J2EiTo;dlu)]+ [ {I + tf E{t-u)dt 



/To 



e+i^(l+t)''(W 



dxdt 



+ lil + t) 



a a \ a 



dxdt. 



Proof, (i) Recall the estimate of u in Proposition 13.11 where fj, is replaced by 6* + 1: 



f I {l + tf+^a{x)u1dxdt < {l + Tof+'^E{To;u)+ f {1 + tf E{t;u)dt 

'To J J To 

+ r hi+tf+^^^^dxdt. 

J To J a (a;) 



We also need Proposition 13.21 with fc = 1: 



/ {l^tf+^E{t-ut)dt < {l + Toy[Y,E{Ta;dlu)]+ f {1 + tf E{t;u)dt 

JTo j^Q JTo 



I To 



i=0 



Since Ai < 1, the claim follows from Proposition 14. If . 
[ii) A simple adaptation of Proposition 13 . II to Ut yields 



/ f{l + tf+^auldxdt < {l + Tof+^E{To;ut)+ f {1 + tf+^E{t;ut)dt 

'To J JTo 

+ r f {i+tf+^^dxdt. 

Jto J a 
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Applying Proposition 13. 21 to the integral of E{t;ut), we obtain 



I fil+tf+^auldxdt < {l+ToY^E{To-dlu)]+ I {1 + E{t;u)dt 



It is clear from Proposition l4.1f M) and Ai, A2 < 1 that aul^ is the main term in the 
upper bound of M'^u. We readily estimate the remaining terms by Proposition 13. 21 
withfc = 2,3. □ 

Remark 4.3. It is straightforward to establish pointwisc estimates of ||a^/^AfM||i2 
and ||a^/^M^u||i2 for large t. Such results will be presented elsewhere. 

5. L°° ESTIMATES IN DIMENSION n = 3 

The most important applications of Propositions 3.4 and 4.2 concern decay 
estimates. We give an example in n = 3, although we can treat all n < 7 due to 
the embedding i/^(R") C L°°(R") for 4 > n/2. It is important to mention that 
the standard Sobolev estimate ^ (I|w||l2 + ||AM||i2) is not suitable, since 

the L°° norm is expected to decay faster than the norm. 

Proposition 5.1. Assume that n = 3 and a, b are -functions. If 

\\{W\nb)u\\L2 < \\b^^^Vu\\L2, 
\\y{a'^^-'u)h2 < \\b'/^Vu\\L2, 

\\^u\\l2 < \\a-'^^divibyu)\\L2, 

for every sufficiently regular u, then 

^ \\a'^^^div{byu)\\L2\\b^/^Vu\\L2. 

Hence ||u||l°° is bounded in terms of E(t;u) and E{t;ut) whenever u is a solution 
of problem {OJ), fO) . 

Proof. We set 

a-^/^div{bVu) = f 

and multiply with u to obtain 

-Au^ = -2 VuP - 2——fu + 2— ■ uVu. 

b b 



The well-known formula 



4:71 J \x — y\ 

which is valid for compactly supported u, yields the estimate 

(5.1) u^ix) < /i--^;7y(^M^)i,,+ /i"^^^ 

J Ky)\x-y\ J \x-y\ 
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It follows from the Cauchy inequality that 

/^^w^r* - 11/11- (/'•(!')'>-'to)i-«'i-'«^H"' 

< ||a->'=div(l)Vo)|U. (J Ma"\y)b-'{y)u{y))fdy\ ' 
where the second factor comes from the Hardy inequality 

y |x - y\-'f{y) dy<ij \Vf{y)\' dy, f E H^K"). 
The assumptions on a and b give 

\V{a'/'{y)b-\y)uiy))\'dy< J b{y)\Vu{y)\Uy, 
so the final estimate becomes 

(5 2) f Ify"^^)', dy<\\a-'/Mivib^u)U4b'^'Vuh.. 

The second integral in (|5.1[) admits similar estimates. From the Cauchy and 
Hardy inequalities, we have 

[ \u(y)\\V\ab-Vu{y)\ , ^ „ ^ 

J \x-y\ 

(5.3) < ||6i/2vy||^,||a-i/2diy(jvM)|U2. 

Adding estimates (|5.2p and (|5.3p completes the proof. □ 



It is easy to find more explicit conditions on a and b which will allow us to use 
Proposition 5.1. This part is about integration by parts and Hardy's inequality. 

Corollary 5.2. Assume that a E and b £ are such that 
(i) a{x) < 1, 1 < b{x), 

(ii) \V{a^^^{x)b-\x))\ + |Vlnfe(x)| < (1 + \x\y\ 
and the matrix with entries 

S' 

(Hi) -^Ab-bj;,^^, i,j^l,...,n 

is non-negative definite. Then the conditions on a and b in Proposition \5.1\ hold; 
hence the L°° estimate also holds. 

Proof. Conditions (i), (ii), and Hardy's inequality readily show that 

\\{Vliib)uh2 < \\b'^'Vuh2, \\V{a'/'b~\)h2 < Wb'^'Vuh. 

for sufficiently regular u. Thus, it remains to verify the condition about ||Au[|^2. 
Let a^^/'^div(6Vw) — f and multiply this equation with Au. Then 

J b{Auy dx = J a^/^fAu dx ~ J (V6 • Vu)Au dx. 

Integrating by parts in the second term on the right side, we obtain 

j (V6 • Vu)Audx ~ 2 y A^I^^P dx — j bx^xjUxiUxj dx. 
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Thus, 



J b{Auf dx = J a^/^fAudx^^ J Ab\Vu\^ dx 

n ^ 

+ ^ / bxix.Ux^Ux. dx. 
* j=i 

We have from assumption {Hi) that 

j b{Aufdx< j a^/^fAudx. 
Now assumption (i) and Cauchy's inequahty imply 

J (Au)' dx< j b{Auf dx < ||/|U.||A^.|U., 
so ||Au||l2 < ||a-i/2div(6Vit)||i2. □ 

6. Application in the homogeneous case c = l,h = 

In this section we will derive explicit decay estimates for (jl.ip when c{x) = 1 
and h{x, t) = 0, i.e. for 

(6.1) utt-d\Y{b{x)Vu) + a{x)ut = 0, x G R", t > 0, 

as an application of the results proven here and in [13] . To state these results we 
introduce the following: 

Hypothesis A. Let a and b satisfy the growth conditions listed in (|2.ip . Then 
there exists a subsolution A{x) which satisfies 

(6.2) div(fe(a;)VA(a;)) > a{x), x G R", 
and has the following properties: 

(6.3) (al) A{x) > for all x, 

(6.4) (a2) A{x) = 0{\x\^-"-'^) for large |a;|, 

a{x)A{x) 

"i'-^'^' b{x)\V A{x)['' 

As we will see below in Theorem 16.11 the quantity fi gives the decay of the 
weighted norm of u, while /i + 1 gives the decay of the energy of u. In several 
cases one can construct explicit subsolutions A which satisfy (al)-(a3), so in these 
cases one can compute the value of /i as given by the formula 

2- a 
^ ~ 2-a-(3' 

These cases are: 

(1) a, b are radial functions. 

(2) a, b are "separable" with the same angular component, i.e. for r — \x\ and 
a; = - € 5"""'^, where S is the unit sphere in R", there exist functions ai, 
&i, and C such that 

a{x) = a(r,cj) = ai(r)C(t^), 
6(a;)=6(r,u;) = 6i(r)C(c^). 

(3) b is radial and o is arbitrary. 



(6.5) (a3) /i:^liminf > 0. 
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(4) a, b are arbitrary (no radial symmetry), a + /3 < 2 and 
yb{x)--^^>boP{l + \x\f-\ 

For a full discussion of these special situations see Section 7 in [13] . 
Theorem 6.1. Assume that a and b satisfy with a, (3 such that 

(6.6) a<l, 0</3<2, 2a + /3 < 2. 

Also, assume that Hypothesis A holds. Then for every 5 > Q the solution of ^6.1\) 
satisfies 

' .S^-S)^a{x)u^ dx < Cs{\\S/uo\\l2 + \\ui\\l2)t^-^', 

e(p-5)^(y2 ^ ^,(^)|vu|2) dx < CsiWVuoWh + \\ui\\l2)t'-^'-^ 

for all t > 1. The constant Cs depends also on R, a, 6, and n. 

An immediate consequence of the above theorem which follows by (al) is 
Corollary 6.2. Under the assumptions of Theorem \6.1[ we have 

(6.7) j a{x)u'' dx < Cs{\\Vu4l2 + ||iii ||i2)^'"^ 

(6.8) E{t-u)= I ul + b{x)\Vu\'' dx<Cs{\\Vu4l2 + \\ul\\l2)t^-^'-\ 



These results in conjunction with the decay estimates presented in Sections 3-5 
yield the following: 

Theorem 6.3. Assume that a, b satisfy i2.1]) and 116.6]) . Then for every S > and 
To > sufficiently large the following inequalities hold: 

(i) Weighted (in time) energy decay: 

(6.9) / {l + t)0+^'^E{t-d^u)dt<{l + TorY.E{To;dlu) + -—^ . 

(ii) Pointwise (in time) energy decay: 
(6.10) 

{l + TorJ2E{n;dlu) + - " 



EiT;d^u)<{l + T) 



+ 5-11 



4=0 

(iii) Weighted kinetic energy decay (the energy associated with the damping): 



(6.11) / auldx < [E{t; u)E{t; Ut)]^'^ 
< (l+T)-"-^ 



1 + nr E ^i^o; diu) + — — -^^ — 
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(6.12) 



(iv) Decay of spatial derivatives (recall that Mu — a ^ [div{bS/ u)]) : 

T 



{l + ty+'a{MuY dxdt 



< 



e+i 



(6.13) 



To 



< 



(1 



Y^E{n-diu) 

.j=o 



rpe + S-fl rpO + S-fJ. 



a{M^uY dxdt 

rpe + S-fJ. 



Y^E{n-diu) 



.i=0 



Remark 6.4. Note from the above decay estimate ()6.10p that the fc-th order energy 
has a polynomial decay of order /i + 1 + 2fc — (5, exactly 2k units below the decay 
of the first order energy as given by (|6.7p . The damping term has the decay rate 
11 + 2 — 5, as it follows from (|6.1ip . which is 2 units below the decay of the 
norm of u given by (|6.8p . Finally, the average decay rates of second and fourth 
order spatial derivatives resemble those of ut and Uu, respectively. Here (1 +t)^^ 
accounts for possibly fast oscillations in b{x)/a{x) as \x\ oo. 

Appendix A. Bounds on the support of solutions 

Recall the non-homogeneous hyperbolic equation with damping 

(A.l) c{x)utt - div(6(a;)VM) + a{x)ut ^ h, x e R", t > 0, 

where the coefficients a, 6, and c satisfy conditions (|2.1[) in the introduction. The 
propagation speed is determined by the ratio c{x)/b{x). We assume the following: 

There exists q € C\ < q'{\x\) < \l x G R". 



(A.2) 



Icjx) 
b{x)- 



A better bound on the propagation speed is given in terms of q satisfying 

|V<z(x)| 



'c(x) 



see Evans "T, but this equation is globally solvable only in special cases. Fortu- 
nately the corresponding inequality is sufficient for propagation speed estimates 
and solvable in most cases. It is also convenient to have bounds depending on |a;|, 
so we look for radial solutions q. 

We consider h G C((0, oo), ^^(R")) and 

(A.3) h{x,t) = if q{\x\) >t + q{R). 

There exist a unique solution u G C((0, oo), H^{W j)nC^{{0, oo), L^{W j) whose 
support is described below. 

Proposition A.l. Let u be a solution of HA.ll) with data {u,Ut)\t=o = (uo,mi), 
such that (-uo,ui) G H^(R'') x ^^(R"). // f^) holds, then uo{x) = ui{x) ^ for 
\x\ > R implies u{x, t) = for q{\x\) > t + q{R)- 

Applying the above result, we can find more explicit bounds when c{x)/b{x) is 
not very small. 
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Corollary A. 2. Let u be a solution of iA . 1\) with data {u,ut)\t=o ~ {uo,ui), such 
that (uo, ui) e iJ^(R") X i^(R"). Assume that b and c satisfy i2.1\) with j3 + ^ < 2. 
There exists q satisfying HA.^) and 

g(|x|) > qo|a;^-(^+^^/^ x G R", 

with qo > 0. Hence, Uo{x) = ui{x) = for \x\ > R implies u{x,t) — for 

't + q{R)V'^'-'-'^^ 



\x\> 



Proof of Provosition HiTTl We can assume that (ito,ui) and h are more regular, 
so M G Ci((0,oo),77i(R")) nC2((0,oo),L2(R")); the general case follows from a 
standard approximation argument. Multiplying (jA.ip with Uf, we obtain 

i {c{x)uj + 5(a;)|V-up)j - div{b{x)utS/u) + a{x)Uf = 0. 

This identity will be integrated over the exterior of a suitable truncated cone. 
Choose q satisfying (jA.2p and consider the conic surface 

M{T)^{{x,t) : T>t>Oandq{\x\)^t + q{R)}, 
which is the lateral boundary of the space-time region 

C(T) = {(a;, t) : T > i > and qi\x\) > t + q{R)}. 
The outer normal to M{T) at (a;,t) is given by the [n -f l)-vector 



We integrate the energy identity over C{T) and apply the divergence theorem: 



{c{x)ul + 6(a;)|Vwp) dx 

q{\x\)>t+q(R) 



a{x)Uf dxdt 



C{T) 



c{x) n bix) ,^ ,o N //I IN X ^\ dS 
Ut+^-^\Vuf + b{x)q{\x\)ut-r-r-'^u^ 



M(T) \ 2 ' 2 ' ' ' " \x\ J yiTMNF 

(c(a;)U( + 6(a;)| Vwp) dx + / hut dxdt, 

qi\x\)>t+qiR) 



JC{T) 

for all T > 0. The integral over M{T) is non-negative, since 

2 2 \x\ 

for all ut and Vu whenever (feg')^ < cb, oi < q' < \fcjb. Thus, 

1 
2 

< \ \ (c(a;)Mi -f 6(x)|Vuop) da; + I hutdxdt. 



(c(x)u^ + 6(x)|VMp) da; + / a{x)u^ dxdt 

qi\x\)>T+q(R) JC{T) 
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Notice that \hut\ < au'f/2 + h'^/{2a). Hence, 
1 



2 



q(\x\)>T+q(B) 



{c{x)u'^ + 6(x)|Vup) dx 



< - / (c(x)u^ + 6(a;)|Vuop) dx + I — dxdt. 

2 Jq(\x\)>q{R) Jc{T) 2a 

The integrands on the right side are due to our assumptions on uq, ui, and h. 
This shows that u{x, T) = if q{\x\) >T + q{R). □ 
Proof of Corollary \ A. SI We already know that u{x,t) = if q{\x\) > t + q{R), 
where (?'(|a;|) < \/c{x)/b{x). It is sufficient to find q, such that 

.'(N) = ^,^, 

where the constants are given in (j2.ip . Solving for q, we obtain 



+ q{Q) = Q, 



6i2-/3-7" 

for some qo > 0. □ 
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